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Abstract 

The peculiar velocities of galaxies distort the clustering pattern anisotropically in 
redshift space. This effect on the statistics of isodensity contours is examined by linear 
theory. The statistics considered in this paper are the three- and two-dimensional genus 
of isodensity contours, the area of isodensity contours, the length of isodensity contours 
in the 2-dimensional slice and the level crossing statistics on the line. We find that all 
these statistics in redshift space as functions of density threshold of contours have the 
same shape as in real space. The redshift space distortion affects only amplitudes of these 
statistics. The three-dimensional genus and the area hardly suffer from the redshift space 
distortion for < Q& -5 / 3 < 1, where b is a linear bias parameter. The other statistics 
are defined in one- or two-dimensional slices of the sample volume and depend on the 
direction of these slices relative to the line of sight. This dependence on direction of these 
statistics provides a way to determine the density parameter of the universe. 

Subject headings: cosmology: theory — galaxies: clustering — galaxies: distances and 
redshifts — gravitation 



*E-mail address: matsu@yayoi.phys.s.u-tokyo. ac.jp 



1. INTRODUCTION 

Redshift surveys of galaxies play essential roles in revealing the structure of our uni- 
verse. If galaxies move purely with the uniform Hubble expansion, redshift surveys would 
tell us the real distribution of galaxies. In reality, however, the peculiar velocities of galax- 
ies distort the distribution in mapping from real space to redshift space. The distortion 
is along the line of sight and the clustering pattern of galaxies in redshift space becomes 
anisotropic. 

There are two characteristic features in the redshift space distortion. On very small 
scales, the random peculiar motions in virializing clusters stretch the shape of clusters 
along the line of sight: the 'finger of God' effect. As a result, the strength of the clustering 
is weaker in redshift space than in real space (e.g., Lilje & Efstathiou 1989; Suto & 
Suginohara 1991; Peacock 1993; Matsubara 1994a). On large scales, the coherent velocity 
field falling in the region with the excess mass make the perturbation enhanced along the 
line of sight, in contrast to the small-scale case (Sargent & Turner 1977; Kaiser 1987; Lilje 
& Efstathiou 1989; McGill 1990). 

The redshift data are to be compared with many theories on the structure formation 
of the universe. The straightforward predictions of these theories, however, are usually 
described in real space. To compare the theories with observations, the various statistical 
measures, such as correlation functions, probability distribution function etc. are used. 
Among such statistics, there is a class of statistics using a smoothed density field which 
cut the noisy property of galaxy distribution. 

For example, Gott, Melott & Dickinson (1986) introduced the topology of isoden- 
sity contours of those smoothed field. The genus G, which is defined by —1/2 times 
the Euler characteristics of two-dimensional surfaces, can be a quantitative measure of 
the topology and was analyzed both in numerical simulations and in redshift surveys of 
galaxies by many authors (Gott, Weinberg & Melott 1987; Weinberg, Gott & Melott 1987; 
Melott, Weinberg & Gott 1988; Gott et al. 1989; Park & Gott 1991; Park, Gott & da 
Costa 1992; Weinberg & Cole 1992; Moore et al. 1992; Vogeley, Park, Geller, Huchra & 
Gott 1994; Rhoads, Gott & Postman 1994). Analytic expressions of the genus for some 
cases are known so far, including Gaussian random field (Doroshkevich 1970; Adler 1981; 
Bardeen et al. 1986; Hamilton, Gott & Weinberg 1986), Rayleigh-Levy random-walk frac- 
tal (Hamilton 1988), union of overlapping balls (Okun 1990) and weakly non-Gaussian 
random fields (Matsubara 1994b). So far these expressions have been derived only for 
isotropic fields. 
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There are other statistics of isodensity contours, which include the 2D genus in two- 
dimensional slices of density field G 2 (Melott et al. 1989), the area of isodensity contours 
iV 3 , the length of isodensity contours in two-dimensional slices of density field N 2 and the 
level crossing statistics Ni (Ryden 1988; Ryden et al. 1989). The analytic expressions of 
all the above statistics for isotropic Gaussian random fields are known. 

The growth of the density fluctuation of the universe on large scales is described by 
linear theory in the gravitational instability picture of the structure formation (e.g., Pee- 
bles 1980). Thus, if the initial fluctuation is a Gaussian random field as is often assumed, 
the statistics of isodensity contours of the density field with large smoothing length should 
obey the random Gaussian prediction. Because the known analytic expression for Gaus- 
sian random fields is for the isotropic field, this Gaussianity test of the initial fluctuation 
should be performed in real space which is not feasible in reality. It is not obvious whether 
or not the redshift space distortion strongly affects statistics of isodensity contours. As 
for the genus, Melott, Weinberg & Gott (1988) found by analysis of iV-body simulations 
that genus is hardly affected by redshift space distortion when the smoothing length is 
larger than the correlation length ~ 5/i _1 Mpc. 

In this paper, the redshift space distortions of statistics of isodensity contours G, G 2} 
iV 3 , N 2} Ni are studied analytically by linear theory of gravitational instability assuming 
that the initial fluctuation is a Gaussian random field. These statistics will be accurately 
determined with the future redshift surveys. Our approach provides the Gaussianity test 
which can be directly performed in redshift space. Moreover, the redshift space distortion 
generally depends on the density parameter of the universe and our formula could in 
principle discriminate this parameter. 

The rest of this paper is organized as follows. In §2, the distant-observer approximation 
to evaluate the redshift space distortion of statistics is introduced. We show directly that 
our definition of this approximation is equivalent to the approximation adopted by Kaiser 
(1987) who gave for the first time the redshift space distortion of two-point statistics by 
linear theory. Then we derive the useful anisotropic statistics in redshift space. The main 
results of this paper, the formula for statistics of isodensity contours in redshift space, are 
derived in §3. We discuss the results in §4. 
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2. FIELD CORRELATIONS IN THE DISTANT-OBSERVER 

APPROXIMATION 



Kaiser (1987) showed that the distortion of power spectrum in redshift space P^(k) 
from that in real space P^ r \h) is given by the simple formula as 

P«(fe)= [l + /p 2 ]V')(*), (2.1) 

where // is the cosine of the angle between the line of sight and the direction of k and 
/(0) = H~ 1 D/D ' 6 , D is the linear growth rate and H is the Hubble parameter. 
The omega- dependence of /(0) is approximately the same in the presence of cosmological 
constant A (see, for detail, Lahav et al. 1991), and we use this approximation, / fi 3 / 5 , 
extensively in this paper. This simplicity of equation (2.1) relies on the approximation that 
the sample volume is distant from the observer. Inhomogeneity of the redshift samples 
closer to the observer is not negligible as well as the anisotropy and prevent to give 
the simple expression as in equation (2.1). When the sample volume is distant from 
the observer, the direction of line of sight is approximately fixed in the sample volume. 
We call this approximation fixing the line of sight as 'distant-observer approximation'. 
Adopting this approximation, the Cartesian coordinates in which the line of sight is fixed is 
convenient for our purpose. The Cartesian coordinates make the calculation of statistics of 
isodensity contours easy. Our distant-observer approximation exactly reproduces Kaiser's 
result (eq.[2.1]) which is derived by first introducing the spherical coordinates and then 
approximating that the sample is distant from the observer. Our approach depends on 
the Cartesian coordinates from the beginning and it would be useful to see directly the 
equivalence of the Kaiser's approximation and our distant-observer approximation. The 
derivation of Kaiser's result in Cartesian coordinates is simpler as we will see in the 
following. 

In the Cartesian coordinates of our distant-observer approximation, we define the 
direction of the line of sight by an unit vector z. Using the line-of-sight component of a 
peculiar velocity field U(r) = v(r) • z, the mapping of the coordinates from the real space 
to the redshift space is given by 

s(r) = r + ^[U(r)-U(0)]. (2.2) 

The observer is placed on the origin of the coordinates, 0. On large scales we are interested 
in, we can relate the number density of galaxies in redshift space and that in real 
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space p(j) by evaluating the Jacobian of the mapping (2.2) resulting in 

p^(s(r))=(l + ±-z.VU(r)y 1 p^(v). (2.3) 



H 

Leaving only linear order in density contrast 6 = p/ p — 1 and peculiar velocity field, this 
relation reduces to 

S (s)(r) = S(;\r)-±2-VU(r). (2.4) 
The peculiar velocity field in linear theory (Peebles 1980) is, in growing mode, 



vlr , 



-tf/VA-^r), (2.5) 



where A -1 is the inverse Laplacian and 8^) is the mass density contrast in real space. In 
the following, 8 g and 8 m are assumed to be proportional to each other. This assumption 
is called linear biasing: 8 g = bS m} where b is the bias parameter which is a constant. The 
relation between the density contrast in redshift space and in real space is, up to linear 
order, 

^(r)=[l + fb-\z-V) 2 A- 1 ]S^(r), (2.6) 

or, in Fourier space, 



1 + fb- 1 



z ■ k\ 2 



k 



? r) (*0, (2-7) 



f2.8) 



which is Kaiser's result. 

In the following, we use the parameters o } and C 3 defined by 

= / ^k«pU{k)W\kR), 
CM = \ \\ dfifi 2 ' (l + /6"V 2 ) 2 , (2-9) 

where V^(x) is the Fourier transform of the window function to smooth the noisy field 
of galaxy distribution. The two popular windows are the Gaussian window Wq{x) = 
exp(— x 2 /2) and the top-hat window Wth^) = 3(sinx — xcosx)/x 3 . We assume that 
the window function is an isotropic function. The rms of density contrast in redshift 
space is given by 

<rW) 2 = C al (2.10) 
We define the following normalized quantities, 

a = Wy A ' = ^P ^ = -^T' (2 - n) 
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(s) 

where 8 R is the smoothed density contrast in redshift space. These quantities obey 
the multivariate Gaussian distribution in linear theory if the primordial fluctuation is a 
random Gaussian field. The multivariate Gaussian distribution is completely determined 
by the correlations of all pairs of variables. For our purpose below, the statistics of 
quantities a, uu (i = 1,2,3; /, J = 1,2) is sufficient. Choosing the coordinates in 
which the line of sight is the third axis, all the correlations among the above quantities 
at some point follows: 
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It is more convenient to consider 



uu = ojij 



a 



(auj u ) 



(2.12) 



instead of u>ij. The new set of variables a, uu distribute as multivariate Gaussian 
and the non- vanishing correlations of these variables are only 
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where 7 = (T^/(a a 2 ). 



3. STATISTICS OF ISODENSITY CONTOURS 

Let us derive the formula of statistics of isodensity contours in redshift space. In the 
following, the primordial fluctuation is assumed to be a random Gaussian field. 

3.1. Genus Statistics 

The genus is minus a half times the Euler number. The Euler number density of 
isodensity contours is evaluated by 

number of maxima + number of minima — number of saddle points (3-1) 

of the contour surfaces with regard to some fixed direction. The expectation value of 
Euler number of the isodensity contours per unit volume is (Doroshkevich 1970; Adler 
1981; Bardeen et al. 1986) 

n^(v) = (S(a - v)6(l3 1 )6(l3 2 )\l3 3 \(u 11 u 22 - ^ 2 2 )) , (3.2) 

where the isodensity contours are defined to be the surface 8$ = va^ s \ This expression 
is valid even for general anisotropic fields. Using new variables ujjj in equation (3.2), the 
following result for the expression of genus in redshift space is derived: 

G^(v) = -\n$\v) = ^^/| (l - |) G"(v), (3.3) 



where, from equation (2.9) 

^ , 1 4- , „ 

5 J 7 



c ° 3 l + ^fb- 1 + -(fb- 1 ) 2 



(3.4) 



and G^ is genus in real space (Doroshkevich 1970; Adler 1980; Bardeen et al. 1986; 
Hamilton et al. 1986) given by 

The redshift space distortion does not alter the shape of genus as a function of density 
threshold and only the amplitude is affected. The dependence of the change in ampli- 
tude is plotted in Fig. 1 (upper panel). The effects of redshift space distortion is small 
for Q& -5 / 3 less than unity. This fact is in agreement with the A^-body analysis of Melott, 
Weinberg & Gott (1988). 
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3.2. 2D Genus Statistics 

The next statistics we consider is 2-dimensional genus. This statistics is defined in 
the 2-dimensional flat plane S in 3-dimensional space. The density field calculated in 
3- dimensional volume defines the high density points in the plane which constitute the 
excursion set on the plane. The 2D genus is defined by the number of contours surrounding 
high density region minus the number of contours surrounding low density regions (Adler 
1980; Coles 1988; Melott et al 1989; Gott et al. 1990). The redshift space is anisotropic 
by the presence of the special direction, line of sight, so the 2D statistics depends on the 
angle 6 s between the plane S and the line of sight. The alternative, equivalent definition 
of 2D genus is useful in the following. For some arbitrarily fixed direction in the 2D 
surface, the maximum and minimum points are defined on the contours. These points 
are classified into upcrossing and downcrossing points with respect to the fixed direction. 
The 2D genus is defined to be 

-(number of upcrossing minima — number of upcrossing maxima 

— number of downcrossing minima + number of downcrossing maxima), 

(3.6) 

of the contour lines with regard to some fixed direction in the plane S. The latter definition 
can be used to obtain the following expression for 2D genus per unit area of the plane: 

GtV, °s) = -\ (% " ^Wi)l& sin #5 + & cos ^K) . (3.7) 

The corresponding expression in case of the isotropic 2-dimensional field is appeared in 
Bond & Efstathiou (1987). From this expression, we obtain 

3 



7 s 




&<'»(,/). (3.8) 



To derive this result, we use u\\ rather than c^n, then regard the variables a, f3 2 sin 6s + 

~ • • • (r) 

/3 3 cos 9s and un as independent variables. The 2D genus in real space G 2 has the 
following form: 

The redshift space distortion again affects only amplitude of 2D genus. The dependence 
on and 6s of the change in amplitude is plotted in Fig. 1 (lower panel). The dependence 
on the direction of the plane 6s for large can be used to determine the cosmological 
parameter by this statistics. 
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3.3. Area of Isodensity Contours 

The area of isodensity contours per unit volume (Ryden 1988; Ryden et al. 1989) is 
given by 

= (s( a - ^A' + Z^ + A?) • (3.10) 
This expression is valid even for general anisotropic fields. Introducing spherical coordi- 
nates for fa, the above expression is calculated to be 

^7 



where is the area in real space given by 



2 




Ci 

Co 




In 





4%) 



(3.11) 



(3.12) 



3 V ; a/37t c 

Again, only the amplitude is affected. The dependence of the amplitude is very weak 
as plotted in Fig. 2 (upper panel). 



3-4- Length of Isodensity Contours in Planes 

As in the case of 2D genus statistics, 2-dimensional flat plane S is considered in 
the statistics of length of isodensity contours. The length of intersections of isodensity 
contours and the plane S was introduced by Ryden (1988). For isotropic density fields, 
this statistics is proportional to the area statistics considered in the previous section. As 
shown below, for anisotropic fields in redshift space, the proportional factor depends on 
the direction of the surface S relative to the line of sight. The angle 6 s between the plane 
S and the line of sight is relevant as in the case of 2D genus statistics. The expectation 
value of length of isodensity contours in the plane S per unit area of the plane is given by 

N^(u, 9 S ) = (S(a - v)^? + (p 2 s\nes + kcose s ) 2 ) . (3.13) 

To evaluate the above equation, note that a, f3i and f3 2 sin #s + /?3 cos 6s are non-correlated, 
independent variables. Introducing the polar coordinates for the latter two variables, the 
following expression is obtained: 

/ f3£i ,\ 2„ \ 



iv 2 (s) Ms; 



1 ) cos 2 9t 



V 



Co V Co 



1 ] cos 2 9 f 



J 



(3.14) 
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where E(k) is the complete elliptical integral of the second kind: 



E(k) = [ V/2 J I - Fsin 2 <j>d<j>, (3.15) 



10 

and is the expectation in real space: 

iV 2 (r V) = I4 r) - (3-16) 

The redshift space distortion affects only amplitude as other statistics considered in this 
paper. The dependence on and 6s of the amplitude is plotted in Fig. 2 (middle panel). 

3.5. Contour Crossings 

Contour crossing statistics is the mean number of intersection of a straight line L and 
the isodensity contours. This statistics of large-scale structure is introduced by Ryden 
(1988) and extensively studied by Ryden et al. (1989) using numerical simulations and 
redshift observations. For isotropic density fields, this statistics is also proportional to 
the area statistics. In redshift space, the density field is anisotropic and this statistics 
depends on the angle 6l between the direction of the line L and the line of sight. The 
mean number of crossings per unit length of the line L is given by 

N[ a \v,Q L ) = (%-i/)|/3 1 sin0 L + # J cos0 L |). (3.17) 

This expression is evaluated by noting that a and fti sin 6l + /?3 cos 6l are non-correlated, 
independent variables, resulting in 



iVfV), (3.18) 



where is the expectation in real space: 



iVfV) = l -Ni r \v). (3.19) 

Again, the redshift space distortion affects only amplitude. The dependence on and 6l 
of the amplitude is plotted in Fig. 2 (lower panel). 
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4. DISCUSSION 



The strength of the effects of redshift space distortion is different according to which 
statistics is focused on. The characteristic point in linear theory is that all the statistics 
in redshift space considered in this paper has the same shapes as in real space as functions 
of density threshold. The redshift space distortion affects only on the amplitude. 

As for genus G(v) and area iV 3 (z/), the redshift distortion of amplitude is small for 
< 06" 5 / 3 < 1 (Fig. 1, 2). This property justifies the comparison of the observational 
redshift data and the theoretical Gaussian prediction in real space [eq. (3.5) and eq. (3.12)] 
at least in linear regime. In the Gaussianity test of primordial fluctuation using genus 
and area statistics, the effect of redshift space distortion can be ignored approximately. 
The similar statistical measure of galaxy clustering, the skewness (6^) / (<5r) 2 induced by 
weakly nonlinear evolution from the Gaussian primordial fluctuation is recently reported 
not to be affected much also by redshift space distortion (Juszkiewicz, Bouchet, k Colombi 
1993; Hivon et al. 1994). 

The direction dependent statistics, 2D genus G 2 , length statistics N 2 and crossing 
statistics Ni are shown to exhibit the dependence on and the direction to define statis- 
tics. In Fig. 3 plotted the direction dependence of these three statistics. The direction 
dependence of amplitude of these statistics is relatively large: for = 1, the am- 

plitude varies more than 20 % while for Q& -5 / 3 = 0, the amplitude does not varies and 
is equal to the one in real space. The direction-dependence depends on Q& -5 / 3 and the 
statistics G 2} N 2 and Ni are three independent indicators to determine the cosmological 
parameters. The redshift space distortion of power spectrum (eq. [2.1]) or two-point cor- 
relation function of Kaiser's result recently used for determining the parameter Q& -5 / 3 
(Hamilton 1992; 1993; Fry k Gaztanaga 1994; Cole, Fisher k Weinberg 1994). Gramann, 
Cen k Gott (1994) introduced the ratio of density gradients ((88^ / dr\\) 2 ) / ((dS^ / dr j_) 2 ) , 
where r_|_, r\\ are spatial component of line of sight and its perpendicular component, as 
a discriminator of Q& -5 / 3 . This ratio is equal to 3C'i/C in linear theory. Our results can 
be used as complementary ways of these observations. 

The purpose of this paper is to give the theoretical relations of statistics of isoden- 
sity contours in redshift space and we did not to apply those results to real data. We 
are planning the applications of our results to data of iV-body simulations and redshift 
observations. 
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FIGURE CAPTIONS 



Figure 1 : Upper panel: amplitude of genus in redshift space relative to that in real 
space. Lower panel: relative amplitude of 2D genus. The angles between the slice 
and the line of sight are S = 0°, 15°, 30°, 45°, 60°, 75°, 90° (from upper line to 
lower line). 

Figure 2 : Upper panel: amplitude of area statistics in redshift space relative to that in 
real space. Middle panel: relative amplitude of length statistics. The angles between 
the slice and the line of sight are S = 0°, 15°, 30°, 45°, 60°, 75°, 90° (from upper 
line to lower line). Lower panel: relative amplitude of level crossing statistics. The 
angles between the crossing line and the line of sight are 6l = 0°, 15°, 30°, 45°, 60°, 
75°, 90° (from upper line to lower line). 

Figure 3 : Relative amplitude of direction-dependent statistics as functions of the angle 
between the line-of-sight and the plane or the line on which the statistics are eval- 
uated. Five cases Q& -5 / 3 = 0.125, 0.25, 0.5, 1.0, 2.0 (from upper line to lower line 
at 6 = 90°) are plotted in each panel. Upper panel: 2D genus. Middle panel: length 
statistics Lower panel: level crossing statistics. 
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